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ABSTRACT 



Aims. In this paper we review some aspects of the theory of magnetic threaded disks. 

Methods. We discuss in particular the equations that determine the position of the inner disk boundary by using different 
prescriptions for the neutron star-accretion disk interaction. We apply the results to several accretion powered X-ray 
pulsars that showed both quasi-periodic oscillations in their X-ray flux and spin-up/spin-down torque reversals. Under 
the hypothesis that the beat-frequency model is applicable to the quasi-periodic oscillations, we show that these sources 
provide an excellent opportunity to test models of the disk-magnetosphere interaction. 

Results. A comparison is carried out between the magnetospheric radius obtained with all the prescriptions used in 
threaded disk models; this shows that none of those prescriptions is able to reproduce the combination of quasi-periodic 
oscillations and torque behaviour observed for different X-ray luminosity levels in the X-ray pulsars in the present 
sample. 

Conclusions. This suggests that the problem of accretion disk threading by stellar magnetic field is still lacking a 
comprehensive solution. We discuss briefly an outline of possible future developments in this field. 
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1. Introduction 

The problem of the interaction between the mag- 
, netic field of a neutron star (NS) and the surround- 
' ing accretion disk has been investigated by a num- 
ber of studi es, with the aid of inagnetohydrodynamic 
simulati ons (iLamb fc Pethickl. 11974 IScharl emamil |1978|; 

Ghosh fc Lambl.ll978HGhosh et al. riiaza iG hosh & Lamb, 

1979a, b: Lovelace et al.'.'1995':'Ro manova et al.l . l2002ll2003l 
1.2004; Ustyugova et al., 2006). Despite several aspects of 
disk-magnetosphere interaction are still poorly understood, 
the idea that the NS magnetic field must penetrate to some 
extend into the accretion disk (due to instabilities leading to 
finite conductivity of the plasma) is now widely accepted. 
This "mag netically threa ded disk mo del" , first develop ed 
in detail bv lGhosh et all ( 1977). Ghos h k Lai^ (|1979allbl 'l 
and later revised bv lWana (,1987, . , 1995) . predicts that, as a 
result of the NS magnetic field threading the disk, a mag- 
netic torque is generated that exchanges angular momen- 
tum between the NS and the disk. The strength of this 
magnetic torque increases in the disk regions closest to the 
NS and exceeds the viscous stresses at the magnetospheric 
radius, i?M, where the disk is terminated. The expression 



//7m-2/7(2GM)-1/7^ 



(1) 
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which holds strictly only for spherically symmetric accre- 
tion, can be considered only a rough approximation when 
the accretion flow is mediated by a disk (in the equation 
above M and /i are the neutron star mass and magnetic mo- 
ment, and M is the mass accretion rate, iLamb fc Pethickl 
Il97l . 

In Sect. [2] we review theories of the threaded disk 
model. Particular attention is given to the differen t mag- 
netic torque prescr iptions of iGhosh fc Lamtj (|1979al lbl) and 
IWand (|1987l 119951) . and the calculation of the magneto- 
spheric radius in the two cases (hereafter GLM and WM, 
respectively). In Sect. [5] these calculations are applied to a 
sample of accretion powered X-ray sources that displayed 
both spin- up /spin- down torque reversa ls and quasi-periodic 
oscillations (QPO. Ivan der Klii . I1995D in their X-ray flux. 
The X-ray luminosity at the onset of the spin- up/spin-down 
transition is used to fix poorly known parameters in the 
equations of the magnetospheric radius in both the GLM 
and WM, whereas the observed QPO frequencies are as- 
sumed to match the predictions of the beat frequency model 
(BFM, see Sect. in order to derive an independent esti- 
mate of Rm- 

We carry out in Sect. [5] a comparison between the mag- 
netospheric radii obtained within the threaded disk models 
and the BFM, and show that neither the GLM nor the WM 
are able to reproduce observations for the whole sample of 
X-ray pulsars considered here. We argue that the problem 
of the threaded accretion disk might still lack a more gen- 
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eral and comprehensive solution, and provide an outline of a 
revision of the GLM that will be presented in a subsequent 
paper (Bozzo et al., 2009, in preparation). 

2. A Review of the Magnetic Threaded Disk Model 

The magnetic threaded disk model was develop ed by 
iGhosh fc LambI (|1979allbh and partly revised by IWand 
(119871119951) 7 under the assumption that the NS is rotating 
about its magnetic axis and that this axis is perpendicular 
to the plane of the disk (the so called "aligned rotator"). 
The model is based on the idea that the stellar magnetic 
field must to some extend penetrate the accretion disk due 
to a variety of effects that prevent this field from being com- 
pletely screened from the disk. Once this occurs, the differ- 
ential motion between the disk, rotating at the Keplerian 
rate fik, and the star, rotating with angular frequency Sis, 
generates a toroidal magnetic field, i?^, from the dipolar 
stellar field component 



BAR) - -v^ 



(2) 



(77 is the screening coefficient, see below). The shear am- 
plification of Bcf, occurs on a time scale rd~|7(ris-i^k)|^^. 
Here n^^iGM/R^y/^ , r2s=27r/P„ is the spin period of 
the NS, h is the disk height, R is the radial distance from 
the NS, and 7~1 parametrizes the steepness of the vertical 
transition between Keplerian motion in the disk and rigid 
corotation with the NS. The finite disk plasma conductiv- 
ity leads to slippage of field lines through the plasma, and 
thus to reconnection of these lines above and below the 
symmetry plane of the disk. This opposes to the shear am- 
plification of the toroidal magnetic field on a time scale 
T^^h/ {S^VAcji)- Here the term £,va4, defines the reconnection 
velocity in terms of the local Alfven speed VA4, and ^~0.01- 
0.1, if the main dissipation effect is the annihilation of the 
poloidal field near the disk midplane, o r if magnetic 
buoyancy is considered. iGhosh fc LambI (|l979aD and lWan^ 
(|l987l[T995[ ) proposed different prescriptions for B^, and in 
the following we discuss their models separately. 

2.1. The Ghosh &. Lamb model 



(3) 



iGhosh fc LambI (|l979a[ ) estimated the toroidal magnetic 
field by equating the amplification and reconnection time- 
scales, i.e. 

Bz iVAz 

In their model, the coupling between the NS and the disk 
occurs in a broad transition zone located between the flow 
inside the NS magnetosphere and the unperturbed disk 
flow. The transition zone comprises two different regions: 

— The boundary layer, that extends from the inner disk 



boundary R^ inward to a distance Ap{<C_R^. In the 
boundary layer the poloidal fleld is twisted at a level 
Bcfj/B^^l, the angular velocity of the disk plasma sig- 
nificantly departs from the Keplerian value, and matter 
leaves the disk plane in the vertical direction and ac- 
cretes onto the star. The currents generated on the disk 
surface by the magnetic field line twisting lead to an 
~80% screening of the NS magnetic field {?]^0.2 out- 
side the boundary layer). 



gi 



— The outer transition zone, where the disk flow is only 
slightly perturbed and the coupling between the disk 
and the star is provided by the residual NS magnetic 
field that survives the screening of the boundary layer. 
Despite this is only a small fraction of the NS magnetic 
field (~20%), the broader extension of the outer transi- 
tion zone makes this coupling appreciable. 

In this model the magnetospheric radius coincides with 
the outer radius of the boundary layer, and is defined as the 
radius at which viscous torques in the disk are balanced by 
magnetic torque s produced by fi eld line twisting within the 
boundary layer. IGhosh fc LambI (fr979a ) found 



— 0-52-Riyr, 



(4) 



where R^ is given by Eq. [TJ 

The total torque N on the star depends on the 
torque iVo=M(GMi?|J)^/^, produced by matter leaving the 
disk at and accreting onto the NS, and the torque 
-^mag=— /tSi BAjBzR^dR generated by twisted magnetic 

field lines threading the disk outside Expressed in an 
adimensional form this torque is 



n = N/No = (iVo+A^magn)/iVo = 1- 



/^gi B^BzR^dR 
M{GMRf^Yn ■ 



(5) 



Beyond i?s, the outer radius of the transition zone, the NS 
magnetic field is completely screened by the disk, and no 
torque is produced. Note that the poloidal magnetic field 
Pz in Eq. [5] differ s from the simpl e dipola r approximation 
of Eq. 12 In fact, iGhosh fc LambI (|l979a| ) calculated this 
component taking into account the effect of the screening 
currents flowing on the disk sur face (see their Eq. 40). As 
a result. [Ghosh fc LambI (|l9791 j) found that n is primarily 
a function of the so called "fastness parameter''^ 



f = njn^{Rt,) = (</i?co) 



3/2 



(6) 



where Rco=1.5xl0^ cm (M/Mq)^/^ p2/3 ^-^e corotation 
radius. For a fixed mass accretion rate, the torque n of Eq.[5] 
can be either positive (the NS spins up) or negative (the NS 
spins down), depending on the NS spin period. In particu- 
lar, for slow rotators (ws<Cl) n~1.4, and the star spins up, 
while increasing u)s, n{u)s) first decreases and then vanishes 
for the critical value Lds=LO^ . For Lds>^^, n becomes nega- 
tive and the star rotation is slowed down by the interaction 
with the accretion disk, whereas for ti;s>0.95 no stationary 
solution exists and steady state accretion is not permitted. 

The critical fastness parameter wf' depends on 
the magnetic pitch angle at the inner disk radius, 
7o=B0(psJ)/B^(i?sJ); IGhosh fc La^ (ll979bD suggested 
7o~l, which corresponds to wf'~0.35 (see their Fig. 4). A 
critical fastness parameter much smaller than unity implies 
the torque on the NS is zero only when the magnetospheric 
radius is well inside the corotation radius and close to the 
compact star (see also Sect.S]). 



^ The explicit functio nal form of n{u 
IGhosh fc LambI (|1979bt ). 



is given by Eq. 7 in 
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2.1.1. The magnetospheric radius in the Ghosh & Lamb 
model 

In the GLM, the magnetospheric radius is given by Eq. U) 
We define the variable Xg\^R^/Rco={ujf)'^^^, and rewrite 
Eq. m in an adimensional form 



2.17;x^/^M-^/^m-i"/2ip-2/3. 



(7) 



Here m is the NS mass in units of IMq. Equation[7]gives 
the ratio Xgi between the magnetospheric and corotation 

radii, for fixed values of ^30, m, Pg, and Mie- Using the 
definition of the critical fastness parameter and defining 



[oj: 



8^)2/3, Eq. [7] translates into 



(8) 



such that the only free parameter is the mass accretion rate 
Mtj-iQ at which a;gi=a;c (i.e. the mass accretion rate at which 
the torque n undergoes a sign reversal). 

In the GLM the magnetospheric radius can be easily es- 
timated from Eq. [51 as a function of the mass accretion rate 
M16, provided the only free parameter Mtrie is somehow 
constrained by observations. This is discussed in Sect.[31 

2.2. The Wang Model 

IWand (|1987f ) suggested that the toroidal field of Eq. [3] is 
overestimated, as the magnetic torque diverges in the limit 
Rs—^oo. Instead of balancing the two time scales and 
Td, he introduced a different prescription for the toroidal 
magnetic field based on Faraday's induction law. Assuming 
that the growth of the toroidal magnetic field is limited by 
reconnection in the disk (see also Sect. I2.1| l. he found that 



± 



£,VAz 



1/2 



(9) 



and proved the amplification of the toroidal field to be 
smaller than previously thought (note that Eq. [5] is the 
square root of Eq. 

In a later studvT iWand ([illl) considered also that the 
growth of the toroidal field might be limited by other 
mechanisms than magnetic reconnection. In case the am- 
plification of the toroidal field is damped by diffusive decay 
due to turbulent mixing w ithin the disk, Tit,—(aflk)~ ^ {a is 
the viscosity parameter of lShakura fc Sunvaevl . ITOTSl here- 
after SS73) and Eq. [5]is replaced by 



B^ ^ j{fls - r^fc) 
Bz ailk 



(10) 



Another possibility is that, in the case of a force- free 
magnetosphere, the winding of the field lines threading the 
disk is limited by magnetic reconnection taking place within 
the magnetosphere itself. In this case 



where 7„aa; < 1- 

At odds with the GLM, the model developed by IWand 
(fl987„ 19.9.5) does not involve the presence of a boundary 
layer: the effect of the screening currents is not taken into 
account self-consistently, and B^ is described by Eq. [21 by 



assuming a constant screening from the inner disk ra- 
dius, i?^, up to the external boundary of the disk (taken to 
be at infinity). Accordingly, the total adimensional torque 
onto the NS is 



n = l- 



/rj^j B^BzR^dR 
M(GMi?J^)i/2 ' 



(12) 



where is given by Eqs.[9l or[Tn]or[TTl and RJ^ is derived 
from the equation 



B^{Rm)Bz{Rm) ^ ] 

M{GMR^y/^ 2{R 



w \3 ■ 

m) 



(13) 



Equation [T31 expresses the balance between the rates at 
which the stellar magnetic field and viscous stresses remove 
angular momentum from the disk. 

The torque given by Eq. [12] is positive for slow rotators 
(wg^l) and negative in the opposite limit, in agreement 
with the results found in the GLM. However, the torque 
vanishes for critical values of the fastness parameter in the 
0. 88-0.95 range, i.e . well ab ove the value a;f'=0.35 predicted 
bv lGhosh fc LambI (|1979aD . In particular, t^^=0.949, 0.875, 
and 0.95 for B^ given by Eqs. [H [TOl and [TTl respectively. 
A similar value is found in appendix [^ where we calculate 
the value of cuc for re gion "C" of a SS73 disk, as opposed 
to region "B" used bv lWand (|l987l '). 

At odds with the GLM, such large values of the crit- 
ical fastness parameters in the WM imply that the mag- 
netospheric radius must lie close to the corotation radius 
when nciO. Therefore, NS spin-down can take place over 
a tiny range of mass accretion rate. This conclusion turns 
out to be nearly independent of the prescription used in 
the WM for the toro idal field. Taking these results into 
account, IWand (|1996f) proposed that a constant magnetic 
pitch at the inner disk boundary, i.e. 70 = B^qB^q, might 
be assumed in Eq. I13i and derived the simplified expression 
for the magnetospheric radius 



R 



M 



(14) 



Here ^~1.357q^^?7''/^, r/ is the screening factor of Eq. [2l 
and R^ is given by Eq. [H 

In Sect. 12.2m we solve Eq.[T31 numerically and compare 
the results with those obtained by using a constant pitch 
angle approximation. 



2.2.1. The magnetospheric radius in the Wang model 

Here we solve Eq. [T^lfor the three different prescriptions of 
the toroidal magnetic field discussed in the previous section. 
We consider first the prescription of Eq. [9l In this case a 
model for the region of the accretion disk that is just outside 
the magnetosphere is required to evaluate the disk h e ight h 
and Alfven velocity waz • In accordance with IWand (Il987t) 
we use the thin disk model of SS73 (see e.g.. iVietril . |2008[ ). 
(11) Using the well known relation h=cl/fti^ connecting the disk 
vertical height h to the sound speed Cs={P / pY^^ (here p 
the matter density and p t he thermal p ressure inside the 
disk), Eq. [3] translates into (jWand . [19871) 



Bz 



7(r!, -f^fc)(47rp)i/2 



Bz 



1/2 



(15) 
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Table 1. Properties of accretion powered X-ray pulsars discussed in the text 



Source 


P. 




fQPO 


iqpo 
(band) 


(band) 


„BFM 
-^QPO 


^QPO 


■^QPO 


Classification^ 




s 


mHz 


mHz 


org s"^ 
(kcV) 


erg s^"^ 
(kcV) 










HerX-li^'-i'^-i^ 


1.24 


806.5 


8 

43 


2.1 X 10^^ 
(0.1 — 200) 


2.1 X 10^^ 
(0.1 — 200) 


0.95-0.98 


0.91-0.97 


0.5 


F 


4U0115-I-63 ".lo-is 


3.62 


276.2 


62 


8 X 10^^ 
(0.1 — 200) 


5 X 10^'^ 
(0.1 — 200) 


0.56 


0.61-0.71 


0.23 


F 


CenX-3 i^'"'!^ 


4.8 


208.3 


35 


1.0 X 10^** 
(0.12 — 100) 


3.3 X 10^^ 
(0.12 — 100) 


0.89 


0.82-0.92 


0.36 


F 


LMCX-4 "'^^ 


13.5 


74.1 


0.65-1.35 


10^» 
(2 — 25) 


2 X 10^** 
(2 — 25) 


0.97-0.98 


0.76-0.88 


0.31 


F 


4U 1626-67 


7.66 


130.6 


40 


1.36x10^'^ 


1.7x10^** 

gQj 


0.83 


0.98-0.99 


0.88 


F/S 


EXO 2030+375 ^'^ 


42 


23.8 


213 


2x10^** 
(1 - 20) 


(1 - 20) 


0.22 


0.86-0.94 


0.4 


S 




42 


23.8 


213 


2x10^* 
(1 - 20) 


2.4 X 10-'"' 
(1-20) 


0.22 


0.43-0.56 


0.14 


s 


A 0535-1-262 ^ 


103 


9.7 


72 


4.3 X lO-''^ 
(20 - 100) 


4.83 X 10^'^ 
(20 - 100) 


0.24 


0.69-0.82 


0.26 


s 


4U 1907+09 


440 


2.3 


55 


6.3 X 10^** 
(2 - 60) 


2 X lO^'' 
(1-15) 


0.12 


0.81-0.92 


0.36 


s 



a: F— Fast rotator, S=Slow rotator. 

References: (1) Parmar ct al. (1989): (2) Angolini ct al. (1989); (3) Finger et al. (1996) : (4) Fritz et al. (2006) : (5) in't Zand et al. (1998) 
(6) Chakrabartv'ct al. (1997); (7) Chakrabartv (1998); (81 Shinoda et al. (1990h o YTamura et al. (1992 ): (10) Campana et al. (2001) 
(11) Woo et al. (1996) : (12) Moon k. Eikenbcrrv (2001a) ; (13) Howe et al. (1983) : (14) Takeshima et al. (1991) ; (15) Dal Fiume et al. (1998) 
(16) Parmar et al. (19991 : (17) Boroson et al. (2000') : (18) Soong fc Swank (1989) ; (191 Burderi et al. (2000') 



Introducing the expressions for the thermal pressure in the 
"B" and "C" regions of the SS73 accretion disk, we obtain 

^-211/80^^ _ ^3/2 ^ 2.72 X 10-V^7-^'7-3"'/'° ■ 

■^,-^''^^y'm'/'pi''''^HIti^ (16) 

and 

^-85/32^^ _ ^3/2 ^ 3^^g^ ^ lO-^v/eFV^"'/'" • 

.M3-o'/Vy4^7/6p85/48^,^63/80^ ^^^^ 

respectively. In the following we refer to these models as 
WMl and WM2, respectively. Here Mig is the mass accre- 
tion rate in units of 10^^ g s~^, ^30 is the magnetic moment 
of the neutron star in units of 10'^'^ G cm"'^, is the mean 
molecular weight, and, in analogy to Sect. 12.1.11 we intro- 
duced the adimensional variable Xw^R-mZ-'^co- With similar 
calculations, we find 

- = 2.38 X 10-3a7-S-V,^o' ' 

■m^/^Pj/^Mie (18) 

and 

^ ~ ■^vi ~ 2.38 X 10 ImaxV ^^30 ' 

•m=^/3pJ/3A/i6, (19) 

by using the prescriptions for the toroidal magnetic field in 
Eqs. (TU] and [Til respectively. In the following we refer to 
these models as WM3 and WM4, respectively. 

Equations [121 [HI [HI ^^nd [TH] give the ratio x between 
the magnctospheric and the corotation radius (we assume 
X < 1), for fixed values of ^, 7, 77, a, /xso, m, Pg, and 
MiQ. Some of these parameters are measured or constrained 



through observations (M, Pg, ™); other parameters are 
still poorly determined by current theory: the values of ^ 
and a (see Sect. 12.2(1 are unce rtain by at least an order of 
magnitude ([King et all . l2007l ) , 77 is in the 0.2-1 range and 
7, Jmax can be larger than 1 (jWand . [l995D . 

In analogy to what we have done in Sect 12.1.11 we use 
here the definition of the critical fastness parameter and 
define <=(w^)2/3. In this case, Eqs. [H [13 [H and [TU 
translate into 

-211/80^1 _ „3/2n1/2 _ 
■^w V-^ "^w ) ~ 

^w -211/80(^ _ 3/2)l/2(^^^/^^^^^^)4/5^ (2O) 

-85/32/-, _ 3/2N1/2 _ 
■^w V-^ -^w / 

^w -85/32(^ - <2'/')^/'(Ml6/Mtrl6)''/''', (21) 

x-'/' - X-' ^ «3-^/2 - <3-2)Mie/Mt,i6, (22) 

and 

- = (x^^-y^ - xr4-2)Mi6/A/tri6, (23) 

respectively. Here x^^=0M6, x^2=0-967, x^3=0.915, 
a;^=0.967, and we defined Mtri6 as the mass accretion rate 
(in unit of 10^^ g s~^) for which Xv,—x'^. In analogy to what 
we found for Eq. ^ Eqs. [101 [H [21 and [Ml show that all 
uncertain parameters cancel out and the magnctospheric 
radius can be easily estimated, as a function of the mass 
accretion rate Miq, provided Mtrw is somehow constrained 
by the observations. This is carried out in Sect. [Hfor the 
sample of X-ray powered pulsars we selected in the present 
study. 

A similar calculation can be applied to the case of 
the constant magnetic pitch approximation. We define 
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0.1 1.0 10.0 100.0 

dM/dt (10" g s ') 



Fig. 1. The ratio x=Rm/ Rco, as a function of the mass 
accretion rate, for different prescriptions of the toroidal 
magnetic fields and Mtri6=3. The solid line is for Eq. [^D] 
(WMl), dotted line for Eq. [2T] (WM2), dashed line for 
Eq. [22] (WM3), and dot-dashed line for Eq. [23] (WM4). 
Thicker lines are the constant magnetic pitch approxima- 
tions corresponding to each of the above prescriptions (note 
that the approximations corresponding to Eqs.[20l[2Tl and 
[23] almost overlap). The triple-dot-dashed line represents 
the GLM (Eq.[SD. 

a^cmp=-RM/F>'co, and rewrite Eq.fMlas 

= A.175^,^^tCM;,'^'M~''/''p-^^\ (24) 

Using the definition of the critical fastness parameter, 
Eq. [23] translates into 

a^cmp = a:r(Ml6/Mt,l6)-'/^ (25) 

where x^=x'^i, a; ^c3^ ^"cii depending on the prescription 
used for the toroidal magnetic field. Note that, in the con- 
stant magnetic pitch approximation, the magnetospheric 
radius does not depend on the prescription used for the 
toroidal field: using a different equation for only af- 
fects the value of the critical fastness parameter that must 
be used in Eq. [21] In Fig. [T] we arbitrarily fixed Mtri6=3 
and compare the values of x—Rm/ Rco, as a function of 
the mass accretion rate, obtained by solving numerically 
Eq. [13] and using the constant pitch approximation for 
all the toroidal magnetic field prescriptions discussed in 
Sect. 12.21 Despite all models approach the same asymp- 
totic behaviour x oc M~^/^ in the limit x<S^l, the constant 
pitch approximation results in a systematic smaller magne- 
tospheric radius, for any considered value of the mass ac- 
cretion rate and prescription of the toroidal magnetic field. 
On the other hand, in the limit x^l, the solution obtained 
with the constant magnetic pitch approximation signifi- 
cantly differs from the numerical solutions obtained with 
Eq. ll3l In particular, the latter results in a magnetospheric 
radius that approaches the corotation radius more gradu- 
ally, while decreasing the mass accretion rate: the range of 
M for which x=Rm/Rco<Xc (and thus the NS spins down), 
is much larger with respect to the range obtained by assum- 
ing a constant magnetic pitch angle (see also Sect. 12. 2p . 
The discrepancy between these results demonstrate that 
the constant magnetic pitch approximation does not pro- 
vide a reliable estimate of the magnetospheric radius in the 



WM. Therefore, we do not use Eq. [^Sin the application to 
X-ray accretion powered pulsars in Sect.[l] For comparison 
we plot in Fig. [1] also Eq. [8] which represents the magne- 
tospheric radius in the GLM. The differences between this 
curve and those derived by using the constant magnetic 
pitch approximation and the WM are due to the different 
values of the critical fastness parameter in the GLM and 
WM (0.35 in the GLM and 0.875-0.95 in the WM; see Eqs.[5] 
and[^. 

3. The beat frequency model 

Besides the threaded disk model, another probe of the 
position of the inner disk radius is offered by observa- 
tions of QPOs in accreting NSs. These timing features 
(jvan der Klisl . 119951 ) have been detected in the X-ray flux 
of a number of astrophysical sources, especially old accret- 
ing NSs and black hole candidates in LMXBs but also in 
young accreting X-ray pulsars in high mass X-ray binaries. 
LMXBs often display a complex variety of simultaneous 
QPO modes, with frequencies ranging from few Hz up to 
'--^l kHz. On the contrary, young X-ray pulsars mostly dis- 
play a single QPO, with a considerably lower frequency 
t'QPO~0.008-0.2 Hz (see Table [1] and references therein). 
Different models have been developed in order to interpret 
the nature of this X-ray variability. The fastest variability, 
manifested through kHz QPOs and timescales of ~ms, must 
be generated by phenomena occurring in the innermost re- 
gions of the accretion disk and reflect the fundamental fre- 
quencies of mo tions in the close v icinity of the compact 
object (see e.g.. lvan der Klid . Il995h. On the contrary, mHz 
QPOs observed in accretion powered X-ray pulsars result 
from variability phenomena occurring farther away from 
the NS (the relevant timescales are of hundreds of seconds). 
In these sources, a magnetic field of order ^10"'^^ G disrupts 
the disk flow at the magnetospheric radius i?M^10^ cm, and 
thus the orbital motion at this radius provides an obvious 
source of variability. However, the involved time scales at 
i?M are a few tens of seconds at the most, and the beat 
between the orbital frequency at this radius and the spin 
frequency of the NS is generally invoked in order to inter- 
pret the observational properties of the slower (mHz) QPOs 

that are observed in these systems. 

According to beat frequency models (BFM. [Alpar et all 
1985; Lamb et al., 1985,), matter from inhomogeneities or- 
biting at the inner disk boundary (Rm) is gradually re- 
moved through the interaction with the neutron star mag- 
netosphere, thus giving giving rise to a modulation in 
the accretion rate and source luminosity. Therefore the 
QPO frequency i^qpo results from the beat between the 
orbital frequency Vovh of the blobs at Rm and the spin 
frequency of the NS i's=27rJ7s, i.e. i'Qpo=t'orb(-RM)-'^s- In 
practice i'orb(-RM) is well approximated with t'K(-RM), *-e- 
the Keplerian frequency at Rm (see below) and the above 
equation can be solved for the magnetospheric radius. This 
gives: 

<,™ = «X10»(^5^)"7^)""™, (26, 

or 

XBFM = 2.2 f (27) 
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where xbfm=^m^'^/-Rco- By assuming that the BFM ap- 
plies, QPOs in X-ray pulsars can be used to probe the physi- 
cal condition of the disk flow at the inner disk radius; in par- 
ticular, Eq. [27] allows for a straightforward estimate of Rm- 
We note that the magnetospheric radius is, by definition, 
the innermost radius at which the disk plasma maintains a 
nearly Keplerian orbit. It might be expected that the beat 
between the disk inhomogeneities (blobs) and the magne- 
tosphere takes place when the former have achieved a sub- 
Keplerian orbital frequency. However, in this case a blob 
would not be centrifugally supported any longer and would 
drift inwards, where the magnetic stresses would rapidly 
bring the blob into corotation with the NS; therefore mod- 
ulated accretion at the beat frequency could not extend 
over several beat cycles, as required to explain the QPO Q- 
fact ors (which range fr om a few to tens in most cases; see 
e.g. Ivan der Klisl . [200? ). Therefore, it can be ruled out that 
disk inhomogeneities giving rise to modulated accretion at 
the beat frequency orbit at substantially slower frequencies 
than Keplerian. On the other hand, a larger orbital fre- 
quency at -Rm than the corresponding Keplerian frequency, 
can be certainly ruled out by the effect of viscosity in the 
accretion disk. 



4. Applications to accretion powered X-ray pulsars 

Here we apply the calculations discussed in the previous 
sections to accretion powered X-ray sources. In particular, 
we selected the sources which displayed QPOs in their X- 
ray flux, as well as evidence for transitions between spin-up 
and spin-down states. For each source, we use the luminos- 
ity measured when QPO were detected (Lqpo) a-nd the lu- 
minosity at which spin-up/spin-down transitions took place 
(Ltr) in order to estimate the magnetospheric radius within 
the BFM (sec Sect. [3]) and magnetically threaded disk mod- 
els (see Sect. [^TTT] and [2.2. ip . respectively. A comparison 
between these estimates of the magnetospheric radius is 
then carried out. In Table [1] we report the relevant values 
of iqpo and Ltr we used, while in appendix |B] we give a 
brief summary of the properties of each source in our sam- 
ple. 

In order to calculate the magnetospheric radius in the 
threaded disk models, we first use the observations of spin- 
up/spin-down transitions. According to the threaded disk 
models, these transitions are the results of changes in the 
sign (from positive to negative) of the torque n acting onto 
the NS. Therefore, the luminosity Ltr can be used to con- 
strain the value of Mtr, at which, according to the models, 
the torque n is expected to undergo a sign reversal (see 
Eqs. [5] and [T^ . In accretion powered X-ray pulsars, the 
conversion between Ltr and Mtr can be obtained by using 
the relation 



Mr36 — 1.3CMt,i6(M/M0)(i?Ns/lO'), 



(28) 



where Ltrse is the X-ray luminosity in unit of 10'^^ erg s^^ 
and <^~1 is an efficiency factor that takes into account, e.g. 
geometrical and bolometric corrections (see later in this 
section) . 

Once Mtrw is determined, the magnetospheric radius 
in both the GLM and WM can be easily estimated as a 
function of the mass accretion rate, since all the uncertain 
parameters cancel out (see Sect. I^TTTT] and [2.2.ip . This is 



shown in Fig. [2] for the X-ray pulsars in our sample (we as- 
sumed a NS mass of m=1.4 and a radius of i?Ns = 10^ cm). 
For each source we plot in the panels of this figure the 
derived values of the magnetospheric radius (units of the 
corotation radius), as a function of the mass accretion rate 
(units of 10^^ g s^^), for the threaded disk models described 
by Eq.EKGLM, triple-dot-dashed line), Eq.dOKWMl, solid 
line), Eq. M (WM2, dotted fine), Eq. ^ (WM3, dashed 
line), and Eq. [23] (WM4, dot-dashed fine). At this point 
we use Eq. [51] and values of Lqpq in Table [1] to calculate 
the mass accretion rate at which QPOs are observed in 
each X-ray pulsar of our sample. The derived mass accre- 
tion rates are represented in panels of Fig. [5] with dotted 
vertical lines. For each source, the intersection between the 
dotted vertical line and the curves representing the mag- 
netic threaded disk models gives the magnetospheric radius 
predicted by these models at the mass accretion rate corre- 
sponding to the Lqpo luminosity. In particular, the inter- 
section with the curves that represent Eq. [201 [H [HI and 
[53] give the range of allowed values of Xvi=Rm/ Rco, *-e. the 
magnetospheric radius in the WM (in unit of the corota- 
tion radius) calculated at the mass accretion rate that cor- 
responds to Lqpq. We indicate this parameter with a^gpo 
in Table [T] Similarly, Xqpq is the value of Xg\—Rf^/Rco for 
which the vertical dotted line intersects the curve from the 
GLM (Eq. [5]) . Finally, we derive for each source the value 
of the magnetospheric radius in the BFM at the mass ac- 
cretion rate that corresponds to Lqpq by using Eq. [^7] We 
call this parameter a;QpQ (a range of values for Xq™ is 
indicated in Table [1] only for those sources which displayed 
more than one QPO frequency). 

With values of a^gpo; ^qpo^ ^^'^ ^qpo hand, the 
GLM and WM can be tested against observations of ac- 
cretion powered X-ray pulsars. Looking at values of these 
parameters in Table [T] we note that the selected sample of 
sources can be roughly divided into two groups. The first 
4 sources (HerX-1, 4U0115-F63, CenX-3, LMCX-4 ) dis- 
played QPO frequencies that, if interpreted in terms of the 
BFM, agree with predictions of the WM. In fact, in these 
cases, 2:qpo ^^'^ -^QPo have similar values, whereas XgpQ 
is typically a factor of 2-3 smaller (for 4U 0115-1-63 and 
LMCX-4 the small discrepancy between a;opo and 



'-QPO 



can be easily accounted for, e.g. by assuming small bolo- 
metric corrections in the X-ray luminosities Lqpq and Ltr, 
see Sect. [U). Values of x close to ~1, as measured for these 
four sources, imply a magnetospheric radius that is very 
close to the corotation radius for the luminosities at which 
QPOs are detected (for example, in the cases of Her X-1 and 
LMC X-4, the magnetospheric and corotation radii differ by 
less than few percent); therefore, in the following we refer 
to these sources as "fast rotators" . Instead, results obtained 
for EXO 2030-^375, A 0535-^262, and 4U 1907-^09 suggest 
the GLM is better suited to account for observations of 
this second group of sources. In this case values of a;QpQ 

are much closer to a^Qpo than Xqpq. However, only for 
A 0535-1-262 a good agreement between the BFM and the 
GLM is obtained. In the other two cases (EXO 2030-^375 
and 4U 1907H-09 ) Xqpq is at least a factor of 2 larger than 

Xqpo {xqpo is a factor of 2-3 larger than a^gpo)- In these 
sources the magnetospheric radius at the mass accretion 
rate corresponding to Lqpq is well inside the corotation 
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Fig. 2. The ratio between the magnetospheric radius and the corotation radius, as a function of the mass accretion 
rate, for sources in Table [TJ The sohd line is for Eq. HD] (WMl), dotted line for Eq. [H] (WM2), dashed hne for Eq. [22] 
(WM3), dot-dashed hne for Eq. [23] (WM4), and triple-dot-dashed line for Eq. E (GLM). In the case of EXO 2030-^375, 
spin- up/spin-down transitions were observed at luminosities in the range ~10^^-2.4x 10'^^ erg s~^, and we represented 
in the figure only those models that provide the higher and lower bound of Rm/ Rco as a function of the mass accretion 
rate. The thinner lines are for Ltr=10^^ erg s"""^, whereas the thicker lines are for Lti— 2.4x lO'^^ erg s~^. In all panels the 
dotted vertical lines represent values of the mass accretion rate corresponding to the luminosity Lqpo at which QPOs 
are observed in each source of our sample. 



radius (x<^l), and thus in the following we refer to them 
as "slow rotators" . The QPO properties of 4U 1626-67 sug- 
gest a magnetospheric radius close to the corotation radius 
(^QPQ— 0.83), like in the case of fast rotators, but they are 
well interpreted within the GLM. This source might thus be 
a sort of "transition object" between fast and slow rotators. 



The conversion in Eq. [25] between observed X-ray lumi- 
nosity and mass accretion rate is affected by several uncer- 
tainties. Besides the NS mass and radius, effects which can 
make C in Eq. [25] differ from unity, such as non-isotropic 
emission and bolometric corrections, should be kept in 
mind. Despite these uncertainties, we note that the results 
derived in this and the next section are virtually insensi- 
tive to variations by a factor of few in Ln- and iqpo • This 
is due to the weak dependence of the magnetospheric ra- 
dius on the mass accretion rate. In all models discussed in 
Sect. [H the steepest dependence of Rm on M is ocM"'^/''; 
therefore, an uncertainty by a factor 2-3 in the X-ray lumi- 
nosity (and thus on M, see Eq. [28]) would cause a 20-30% 
change in the magnetospheric radius at the most. 



4.1. Two case studies: 4U 1626-67 and CenX-3 

In this section we discuss further the cases of 4U 1626-67 
and CenX-3, for which detailed studies of the long term 
variations of the QPO frequency with the X-ray flux were 
recently published. 

In particular, iKaur et all (|2008D studied QPOs in 
4U 1626-67 at different X-ray fluxes and compared the ob- 
served frequencies with those calculated by using the BFM 
and with the magnetospheric radius determined based on 
the GLM. These authors noted a discrepancy between the 
observations and the predictions, and argued that the BFM 
might not apply to this source. In Fig [3] (upper panel) we 
show the same calculation, but included also the QPO fre- 
quencies estimated by using the BFM with the magneto- 
spheric radius as determined in the WM. In this plot the 
solid line corresponds to the QPO frequency predicted by 
the BFM when Eq.l^fWMl) is assumed for the magneto- 
spheric radius, the dashed line is for Eq. [2I](WM2), dotted 
line for Eq. [22] (WM3), and the dash-dotted hne for Eq.[23] 
(WM4). The triple-dot-dashed line represents the QPO fre- 
quencies predicted by the BFM when the magnetospheric 
radius is calculated ac cording to the GL M ( Eg. [8]). We se- 
lected those data from lKaur et al.l (|2008f ) and lKrauss et al.l 
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Fig. 3. Predicted QPO frequencies in the BFM when the 
magnetospheric radius is calculated according to Eq. [^D] 
(WMl, solid line), Eq. [2T] (WM2, dashed line) Eq. ^ 
(WM3, dotted hne), Eq. [22] (WM4, dash-dotted line), and 
Eq.[5](GLM, triple-dot-dashed line). The upper panel is for 
4U 1626-67, whereas the lower panel is for CenX-3. In both 
cases QPO frequencies inferred from observations of these 
sources at different X-ray fluxes are shown to gether with 
measu r ement errors. These data taken from Kaur et ahl 
(jlOOi), iKrauss etal] (|2007D . and iRaichur fc Paul llOOl); 
we used Eq. [2H1 to convert between X-ray fluxes and mass 
accretion rates. 



(j2007t ) for which QPO frequencies and X-ray fluxes were 
measured simultaneously, and used Eq. [25]to convert these 
fluxes into mass accretion rates. In the upper panel of Fig. [3] 
we show that, due to the different prescriptions available for 
the magnetospheric radius (i.e. the GLM and the WM), 
the region of predicted QPO frequencies in the BFM, as a 
function of the mass accretion rate, is very broad and all 
observational measurements lie within this region. 

The lower panel of Fig. [3l show the ca se of CenX-3. 
We used data from iRaichur fc Paull ([lOO^)- In this work, 
the authors showed that the QPO frequency of this source 
has almost no dependence on the X-ray flux. By using the 
GLM to calculate the magnetospheric radius, they argued 
that, if the BFM applies, then the long term X-ray inten- 
sity variations of CenX-3 are likely due to obscuration by 
an aperiodically precessing warped disk, rather than being 
related to changes in the mass accretion rate (and thus lo- 
cation of the inner disk radius). In fact, in the latter case 
the QPO frequency would be expected to vary according 



to Eqs. [27l and l4l However, our calculations show that all 
measured QPO frequencies lie inside the region spanned by 
different magnetically threaded disk models. We conclude 
that the observations of 4U 1626-67 and CenX-3 do not 
support simple applications of either the GLM or WM to 
the BFM. We further comment on this in Sect. [S] 



5. Discussion 

The results obtained in the previous section indicate nei- 
ther the GLM nor the WM, when used in conjunction with 
the BFM, are able to reproduce the range of observations 
discussed here for the entire sample of X-ray pulsars. We 
also note that for all sources in Fig. [21 the magnetospheric 
radius in the GLM turns out to be somewhat smaller than 
that de rived by usin g the WM. This point was discussed 
also by IWan j (jl996l ). who suggested that the reason for 
this disagreement resides in the different prescription of the 
toroidal field used in the two models: the assumed in 
the GLM implies a larger magnetic torque that spins down 
the NS more efficiently and reduces the value of the critical 
fastness parameter. As a consequence, the GLM magneto- 
spheric radius is located closer to the NS (see Fig. [2]). We 
showed here that the magnetospheric radius predicted by 
the GLM is still too large to account for observations of 
the QPOs over the entire sample of the slow rotators (see 
Table [IJ . This remains true even when the GLM is revised 
to include a more accurate prescription of the toroidal mag- 
netic field, whic h leads to larger values of the critical fast - 
ness parameter (IWand . [19871 iGhosh fc La^ . 11991119951) . 

At odds with the GLM, the magnetospheric radius in 
the WM approaches the corotation radius more gradually 
as the mass accretion rate decreases, a result that seems 
to account for observations of fast rotating sources (see 
Table [J). However, in the cases of 4U 1626-67 and CenX-3, 
for which detailed studies of the long term variations of the 
QPO frequency with the X-ray flux are available, the WM 
is not able to reproduce the observations. It was also noted 
that the treatment of the NS poloidal field screening by 
currents fiowing o nto the disk sur f ace in this model might 
be oversimplified fG hosh fc La^ . Il992l [19951) . 

Furthermore, in Sect. 12.21 we pointed out that an im- 
portant caveat in the WM is that the interaction between 
the accretion disk and the NS magnetic field takes place 
in a similar fashion over the whole accretion disk. This is 
at odds with the GLM that predicts the strong coupling 
between the NS and the disk takes place mostly within a 
small boundary layer, such that this region alone deter- 
mines the position where the disk terminates (i.e. Ry\). On 
the other hand, the theory of the boundary layer envisaged 
in the GLM might not be applicable to fast rotators, being 
the radial extent of the boundary layer of the same or- 
der of the separation between the magnetospheric and the 
corotation radii in these casefl Some works have investi- 
gated the i mportance of the boundary la yer in the threaded 
disk model (iLi fc Wand.ll996HLi fc Wic kramasinghe. 199^ 
iLi fc Wand . 119991 : lErkut fc Alparl . [2004R ,Li fc WaQ 



^ Note that I Ghosh fc LambI l|l979bl ) showed steady state ac- 
cretion might not be allowed in their model for ojc^O.OS. 

^ "Torqueless accre tion" l|Li et al.l . Il996l : 

ILi fc Wickramasinghd . 119971 1 was not considered here since that 
mechanism is unlikely to b e applicable to accretio n powered 
X-ray pulsars (jWand . lT997bl : iRomanova et al.l . l2002l 'l. 
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(|1996f ) suggeste d that t here e xists an uncertainty by a fac 
tor of 4 in the IWand ()1987f ) equation defining th e ma; 
netospheric radius; using this result, iLi fc Waii3 ( |199' 
demonstrated that the boundary layer in NS accreting bi- 
naries can survive the destructing action of the NS mag- 
netic field down to a radius ~0.8i?M- The boundary layer 
might thus be si g nific antly larger than previously thought 
([Ghosh fc Lambl . Il979a ). However, the derived corrected 
value of the critical fastness parameter (^0.71 - 0.95) does 
not differ much from previous estimates (|Wangl . ll99,5 r> and 
the problem of slow rotatin g sources remain s open . Similar 
results were obtained by lErkut fc Alparl (|2004f ). which 
demonstrated that the width of the boundary layer might 
be a strong function of the fastness parameter: they found 
that broad "boundary layers" are expected for spinning-up 
sources, whereas much reduced boundary layers should be 
expected for sources in a spin-down state (these boundary 
layers are typically a factor ~6-60 less wider than those 
found for spinning-up sources). However, a general analyti- 
cal equation for the magnetospheric radius cannot be easily 
derived, due to the presence of few additional parameters 
in their model. Br oad boundary layers we re found also in 
the simulations bv iRomanova et ail (|2002l ). These authors 
found a reasonable agreement with the predictions of the 
GLM, with the inner region of the disk behaving like a 
boundary layer, while the outer region is only partially cou- 
pled with the magnetic field of the star. These numerical 
simulations suggested a critical fastness parameter of ~0.6. 
Despite this value is in between the values obtained within 
the GLM and the WM, it cannot account for observations 
of both fast and slow rotators. Our results in Table [T] imply 
that the critical fastness parameter cannot be constant for 
all these sources. A more general solution for the magnetic 
threaded disk model might be found in the future in which 
the WM and the GLM give the limiting cases of fast and 
slow rotation, respectively. 

The present study suggests that all the discussed lim- 
itations of both the WM and GLM might be the reason 
why none of these models is able to reproduce the com- 
bination of QPO and torque behaviour observed at differ- 
ent X-ray luminosity levels in the X-ray pulsars considered. 
Alternatively, the BFM might not be applicable to (all) 
QPOs observed from X-ray pulsars. 



6. Conclusions 

We showed that, if the BFM applies to the QPOs of X-ray 
pulsars, then the GLM and WM cannot completely account 
for observations of the sources in our sample. Instead, tak- 
ing into accounts results in Table [TJ we noted that these 
sources can be divided into two classes: 

— Fast rotators, for which the Keplerian velocity of matter 
at the inner disk radius, as inferred from the application 
of the BFM to the observed QPO frequency, is close to 
the rotational velocity of the star. In this case, the mag- 
netospheric radius inferred from the WM and the BFM 
predict QPO frequencies which seem in good agreement 
with the observations. However, we showed in Sect 14.11 
that at least in the cases of 4U 1626-67 and CcnX-3, 
the WM is not able to reproduce observations of QPO 
frequencies at different X-ray fluxes. 

— Slow rotators, for which the Keplerian velocity of mat- 
ter at the inner disk radius, as inferred from the appli- 



cation of the BFM to the observed QPO frequency, is 
well above the spin velocity of the NS. In this case, the 
magnetospheric radius derived from the BFM is less dis- 
cordant with the predictions of the GLM. In fact, only 
for A 0535-1-262 a good agreement between the GLM 
and the observations is obtained. For slow rotators the 
WM give a magnetospheric radius that is at least >2-8 
times larger than that derived from the BFM. 

We conclude that either a more advanced theory of mag- 
netically threaded disks is required, or that the BFM does 
not apply to (all) QPOs observed from X-ray pulsars. 

Appendix A: Calculation of the torque for region 
"C" of the SS73 accretion disk model 



In order to evaluate p and h in Eas.[giand [T^IWan3 (|1987l ) 
considered only the "B" region of the SS73 accretion disk, 
i.e. the gas- pressure dominated region where electron scat- 
tering gives the main contribution to the opacity. Here we 
carry out the same calculation by using the expressions of 
p and h that are appropriate to region "C" of the SS73 
accretion disk model (where the main contribution to the 
opacity is pr ovided by free- free absorption). According to 
IVietril (|2008f) . the thermal pressure of disk matter has a ra- 
dial dependence p = _pm(-Rm/-R)^"^^^, where the subscript 
M denotes quantities evaluated at the inner disk radius. 
From Eq. [H] we get 
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where the same notation as that of Sect. [2] is used. Setting 
y = {R/Rcof/'^ and using Ea. lA.2l into Eq.[n]we obtain 
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By numerically evaluating integrals in the above equation 
we find LOr — 0.95. 



Appendix B: Values of Ltr and -Lqpo for accretion 
powered X-ray sources 

Here we briefly summarize the relevant observations of the 
accretion powered X-ray pulsars considered in Table (TJ in 
order to explain values used for the luminosities Ltr and 

-^QPO- 

— HerX-1: HerX-1 is one of the best studied X-ray 
binary system. It consists of a ~1.24 s spinning 
NS and a A/F companion (the orbital period is 
1.7 day). The X-ray flux displays a regular modu- 
lation at a 35 day period, that has been associated 
to the precession of a highly warp ed accretion disk 
that periodically obscures the NS (iPett erson L 11975 : 
Choi et all. 119941: IWilson et al.1. Il997t iDa l Fiumc et al. , 
19981 : iParmar et all . 119991 : iKlochkov eTal., . .2007,1 . This 
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suggested that transitions between high ("main-on 
states" ) and low luminosity states (the "anomalous" low 
states) of this source are to be interpreted as due to lo- 
cal obscuration phenomena, rather than large changes in 
the mass accretion rate. Evidence in favor of this inter- 
pretation has been recently obtained th r ough detailed 
phase-resolved spectroscopy (jZane et al.l . l2004l ) as well 
as observations of X-r ay heating of the companion star 
(jBoroson et all |2000( ). Since QPOs and spin- up/spin- 
down transitions were observed during both high and 
low luminosity states, for the purpose of this paper 
we assume Ltr=iQPo=2.1 x 10^^ erg s~^ (a distance of 
5 kpc is considered), where the latter is the typical main- 
on state luminosity. 

4U 0115+63: 4U 0115-1-63 is a binary system hosting a 
3.6 s spinning NS orbiting a Be companion (the or- 
bital per iod is 24.3 day). The d i stance to the source is 
~8 kpc (jNegueruela &: Okazakil . l200lh . An HEAO ob- 
servation caught this source in outburst (the typical out- 
burst luminosity is ~8xl0^^ erg s~^. ICampana et al.l . 
[2001), and a prominent peak in the power spectrum 
of the X-ray lig ht curve was detected at 62 mHz 
(|Soong k Swankl . ll989D . iTamura et all (|1992[ ) reported 
that the spin-up trend the source usually displayed 
while in outbursts reversed during lower luminosity 
states. These typicall y occurred at SxlO^^^ erg s^^ 
(jCampana et al.lboOlh . 

CenX-3: CenX-3 is a high mass X-ray binary with a 
spin period of ~4.8 s and an orbital period of '--^2.1 day. 
The companion star is an 0-typc supcrgi ant and the dis- 
tance is estimated to be ~8 kpc (Bur deri et al.l . l2000l 
and references therein). The presence of a ~35 mHz 
QPO in the power spectrum of this source was first 
discovered bv iTakeshima et all (jl991h . after the source 
egress from an X-ray eclipse. The typical X-ray lumi- 
nosity was determined with BeppoSAX, and is of order 
~1.0xl038 erg s-i ( 0.12-100 keV). CenX -3 has a secu- 
lar spin-down trend ijBildsten et al.l . [T997[ ). but episodes 
of spin reversal were found to occur at a luminosity that 
is typically a factor of ;^3 below that of t he post-eclipse 
high luminosity state (jHowe et al.Ul983[ ). We thus con- 
sidered in Table [T] that Lqpo is equal to the typical X- 
ray luminosity observed in the post-eclipse egress state, 
whereas Ltr'^l/3LQPo. 

LMCX-4: LMCX-4 is an accretion powered X-ray pul- 
sar with a spin period of ~13.5 s and an orbital pe- 
riod of ^1.4 day. The X-ray luminosity of the sys- 
tem varies with a periodicity of ~30.3 day, alternat- 
ing between high (~2xl0'^* erg s~^) and low (a fac- 
tor of ~60 below) lu minosity sta t es (fo r an estimated 
distance of ~50 kpc, IWoo et all . |1996[ ). This period- 
icity have been attribute d to the effect of an obscur- 
ing tilted accretion disk (jMoon fc Eikenberr^ l200lU 
and references therein). During the high states, spin 
torque reversals were repeatedly observed, whereas dur- 
ing episodes of very bright flares (~10^^ erg s~^, 2-25 
kcV) QPOs were detected at frequencies in the ~0.65- 
1.35 mHz range (these w ere inte rpreted within the BFM 
by Moon fc Eik cnberr^, l2001bl ) . 

4 U 1626-67: 4U 1626-67 is a low mass X-ray binary host- 
ing a ~7.7 s spinning neutron star in a ^^42 min orbit 
around a ^^0.004 M0 companion star. QPOs were de- 
tected in the X-ray observa tions of this source more than 
once (for a review see, e.g. iKaur et al.L[2"00a) . A torque 



reversal was observed bv lChakrabartv et al.l (|1997f ). who 
also estimated a source distance of •^3 kpc. For a review 
of the flux history o f 4U 1626-67, we reefer the reader to 
iKrauss et ahl (|2007D . 

- EXO 2030+375: EXO 2030-1-375 is a Be X-ray tran- 
sient with an orbital period of ~46 day, hosting a 
~42 s spinning n eutron star, and loc ated at a dis- 
tance of 7.1 kpc (jWilson et al.l . l2002f) . In 1985 this 
source underwent a bright outburst (peak luminosity 
~2xl0 ^^ erg s~M, arid a Q PO at ~213 mHz was de- 
tected (jAngelini et al.l . [T989I ) . Spin- up/spin-down tran- 
sitions were observed more than once, at luminosities in 
the -10^8-2.4x1036 erg s'^ (1-20 keV) range. 

- A0535+262: A0535-h262 is a ~103 s X-ray pulsar, or- 
biting a 09.7 companion star (the orbital period is 
'^lll day). QPOs and spin reversals were best observe d 
during the giant outburst in 1994 (Fing er et al.l . [T996| ). 
This outburst was detected with BATSE in the energy 
20-100 keV range, and the flux at the peak of the out- 
burst was ~6 Crab. Observations of this outburst at 
lower energies (<20 keV) were no t availa ble, but, based 
on previous results, iFinger et al.l ()1996[ ) estimated the 
2-10 keV flux might not be larger than 2 Crab. Due to 
the uncertainties in this estimate we have not corrected 
values reported in Table[T]for this factor. As discussed in 
Sect. Uan uncertainty of ^^30% on the luminosity used 
to derive the position of the magnetospheric radius can- 
no t affect much our re sults. Note that the distance used 
bv lFinger et al.l (|1996[ ) to con vert the observed fl ux into 
an X-ray luminosity is 2 kpc ([Steele et al.l . [1998') . 

- 4U 1907+09: 4U1907-f09 is a ~440 s spinning NS 
in a ~8 day orbit around a supergiant companion. 
QPOs were discovered during an hour long flare at 
~6. 3x10 ^6 erg s~^ (2-60 keV, and an assumed d istance 
of 5 kpc, lin't Zand etlll . 119981: ICox et _a l.'. 2005) . After 
the d iscovery of the pulsations bv [Makishima et al.l 
([1992f) . the NS in 4U 1907-^09 exhibited a steady spin- 
down for about 20 yrs. This trend changed in 2006, when 
a torqu e reversal was obs erved at ~2xl0'^6 erg s^^ (1- 
15 keV, iFritz et al.l . [20061) . 
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